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SPECIFICATION 

MAGNETIC RESONANCE IMAGING APPARATUS 

Technical Field 

The present invention relates to a magnetic resonance 
imaging apparatus (hereinafter referred to as "MRI apparatus") 
for imaging a tissue structure inside an object to be examined 
by utilizing a magnetic resonance phenomenon, more 
particularly to an MRI apparatus suitable for parallel imaging for 
rapidly imaging the object using a plurality of receiving coils. 

Background of the Invention 

In an MRI apparatus, a pulse sequence is sequentially 
executed while a phase encoding amount given to a nuclear 
magnetic resonance signal (hereinafter referred to as "NMR 
signal") is varied, and an echo signal necessary for 
reconstructing one image is thus acquired. Therefore, imaging 
time depends on time of repetition (TR) and the repetition 
number of times of pulse sequence. When the rapid imaging is 
performed, a multi-echo pulse sequence in which plural echo 
signals can be measured by executing the pulse sequence one 
time and a pulse sequence in which the TR is shortened to 
several ms to several ten ms are generally used. 

Meanwhile, when a region around a heart, such as a 
coronary artery, is imaged, it is necessary to perform imaging 
more rapidly than in the imaging based on the above pulse 
sequence. To realize such rapid imaging, a rapid imaging 
method called "parallel imaging method" is proposed. 



The parallel imaging method is a method of reducing the 
repetition number of pulse sequence by executing the pulse 
sequence using a plurality of receiving coils while phase 
encoding steps are thinned out at even intervals, and thus 
shortening the imaging time. When a measured signal 
obtained by performing the parallel imaging pulse sequence is 
simply image-reconstructed, an aliasing artifact is generated in 
the reconstructed image. The aliasing artifact is generated 
due to a difference between a field of view (FOV) of image to be 
obtained and that in which signals are actually measured in the 
parallel imaging, and it can be removed by performing a matrix 
calculation using a pre-calculated sensitivity distribution of 
each receiving coil. 

Generally, in the parallel imaging method, the imaging 
time can be shortened depending on the number of receiving 
coils. That is, it is possible in principle to thin out the number 
of phase encoding steps matrix as much as the number of 
receiving coils from the number of phase encoding steps in the 
usual imaging. As a result, for example, the imaging time can 
be shortened 1/N times, where the number of receiving coil is N, 
in comparison with the imaging performed with the usual phase 
encoding step number. 

However, in the parallel imaging method, a sensitivity 
distribution of each receiving coil is needed for removing the 
aliasing artifact generated in the image. To obtain the 
sensitivity distribution of each receiving coil, it is necessary 
before imaging to execute a pulse sequence for acquiring a 
sensitivity distribution of receiving coil in the state where the 
object is put in an imaging space and wearing the receiving coil, 



in addition to the pulse sequence for obtaining an examination 
image of the object in which the phase encoding steps- are 
thinned out. Moreover, when imaging is performed on many 
slices on various positions in an imaging portion of the object to 
conduct diagnosis, it is necessary to execute the pulse 
sequence for the sensitivity distribution acquisition on every 
slice position to be imaged. Therefore, there is a problem that 
the total imaging time in the parallel imaging becomes very 
long. 

Techniques of shortening the total imaging time in the 
parallel imaging includes one in which the sensitivity 
distribution of receiving coil is acquired by performing 
three-dimensional measurement. However, because the 
sensitivity distribution data are similarly acquired on every 
imaging slice position, the effect of shortening the time is small, 
and it is desired that the time be further shortened. 

Furthermore, when the imaging is performed again on a 
different slice position after measuring the sensitivity 
distribution of receiving coils in advance of picturing and 
imaging the object, it is necessary to re-measure the sensitivity 
distribution of receiving coil on the present imaging slice 
position. 

A first object of the present invention is to provide an 
MRI apparatus which can shorten the total imaging time when 
plural slices of the imaging portion of the object are pictured 
using the parallel imaging method. 

A second object of the present invention is to provide an 
MRI apparatus which can acquire the sensitivity distribution of 
receiving coil necessary for removing the aliasing artifact within 



a short time when the object is imaged using the parallel 
imaging method. 

A third object of the present invention is to provide an 
MRI apparatus which can acquire the sensitivity distribution of 
5 receiving coil necessary for removing the aliasing artifact in a 
simple manner when the object is imaged using the parallel 
imaging method. 

Summary of the Invention 

10 To achieve the above objects, an MRI apparatus 

according to the present invention includes a plurality of 
receiving coils for receiving a magnetic resonance signal 
generated from an object to be examined, sensitivity image data 
acquiring means for acquiring sensitivity image data of each of 

15 the plural receiving coils from plural slice positions separated 
from each other on the object by executing a first pulse 
sequence, means for acquiring examination image data of each 
receiving coil from the plural slice positions adjoining each 
other on the object by executing a second pulse sequence using 

20 the plural receiving coils in a which phase encoding matrix in a 
k space is thinned out, and an artifact removing means for 
generating sensitivity distribution data of the plural receiving 
coils on slice positions of the examination image data on the 
basis of the plural sensitivity image data and removing an 

25 aliasing artifact in the examination image using thus generated 
sensitivity distribution data of receiving coil. 

In the image data of a portion appearing as the aliasing 
artifact in the image, a data component of one portion is 
wrapped and superposed on that of another portion. 
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Accordingly, by forming a simultaneous equation in accordance 
with a correlation among the image data of the plural receiving 
coils, the sensitivity distribution data of each receiving coil, and 
the data components on plural points of the imaging cross 
5 section and solving these simultaneous equations on the basis 
of matrix calculation, the superposing data component at the 
wrapped portion is separated. The aliasing artifact on the 
image is removed by expanding thus separated image data. 

In this manner, to remove the component of aliasing 

10 artifact in the image data acquired by executing the second 
pulse sequence in which the phase encoding steps in the 
measurement space, i.e. the k space are thinned out, the 
sensitivity distribution data on each slice position of image data 
are necessary. These sensitivity distribution data are 

15 calculated from the sensitivity image data acquired by executing 
the first pulse sequence for acquiring measured data of central 
portion (low-frequency portion) in the k space, whereby the time 
for acquiring the sensitivity distribution of receiving coil can be 
shortened. 

20 In the MRI apparatus according to the present invention, 

the slice number of the sensitivity image data acquired in the 

first pulse sequence is n in number, which is smaller than m in 

i 

number of examination image data including the aliasing artifact 
component acquired in the second pulse sequence. When 
25 sensitivity distribution data on a slice position do not exist in 
the image data for m slices, they are generated by performing 
interpolation using the above acquired sensitivity distribution 
data. Therefore, because the repetition number of the first 
pulse sequence for acquiring the sensitivity distribution of 
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receiving coil can be reduced, the total time of repeatedly 
executing the first pulse sequence is shortened, and the total 
imaging time in the parallel imaging including the first pulse 
sequence and the second pulse sequence can be thus 
5 shortened. 

Further, the MRI apparatus according to the present 
invention includes means for acquiring sensitivity image data of 
each receiving coil for plural slices on slice positions different 
from each other by executing the first sequence and acquiring 

10 aliasing image data of each receiving coil for plural images by 
executing the second sequence, and means for generating the 
sensitivity distribution data of each receiving coil on the slice; 
positions of the above sensitivity image data on the basis 
thereof and removing the aliasing artifact generated in the 

15 examination image on the basis of the sensitivity distribution 
data on the slice positions of the aliasing image data. 

Further, in the MRI apparatus according to the present 
invention, the plurality of receiving coils may include a first 
receiving coil and a second receiving coil. In this case, it is 

20 preferable that the sensitivity distribution of the first receiving 
coil is more uniform than that of the second receiving coil. 
Thus, by using as a reference an image having a uniform 
sensitivity distribution obtained by the first receiving coil, the 
sensitivity distribution of the second receiving coil can be 

25 accurately calculated. 

Further, in a structure where the first receiving coil does 
not have a uniform sensitivity distribution, an image having a 
more uniform sensitivity distribution than that generated on the 
basis of the measured data acquired by a single receiving coil 



is generated by combining the measured data obtained by using 
a plurality of receiving coils. By using as a reference* this 
image having a uniform sensitivity distribution, the sensitivity 
distribution of each receiving coil can be accurately calculated. 

Brief Description of the Drawings 

Fig . 1 is a diagram showing a procedure ,of signal 
processing of a signal processing unit in an MRI apparatus 
according to a first embodiment of the present invention. Fig.2 
is a diagram showing a schematic structure of the MRI 
apparatus according to the first . embodiment of the present 
invention. Fig. 3 is a diagram illustrating a receiving coil in the 
MRI apparatus according to the first embodiment of the present 
invention. Fig. 4 is a diagram showing a schematic structure of 
the multiple RF receiving coils and a reception unit in the MRI 
apparatus according to the first embodiment of the present 
invention. Fig. 5 is a chart illustrating a pulse sequence used 
in the MRI apparatus according to the present invention. Fig. 6 
is a diagram illustrating an aliasing artifact to be removed in the 
parallel imaging method in the MRI apparatus according to the 
present invention. Fig. 7 is a diagram illustrating image data 
acquired in the pulse sequence and those generated by a slice 
interpolation processing in the MRI apparatus according to the 
first embodiment of the present invention. Fig. 8 is a diagram 
showing a procedure of signal processing in the signal 
processing unit in the MRI apparatus according to a second 
embodiment of the present invention. Fig. 9 is a diagram 
illustrating sensitivity distribution data generated by a 
sensitivity distribution calculation processing and a slice 
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interpolation processing and image data acquired in the pulse 
sequence in the MRI apparatus according to the second 
embodiment of the present invention. Fig. 10 is a diagram 
illustrating a receiving coil in the MRI apparatus according to a 
third embodiment of the present invention. Fig. 11 is a diagram 
showing a schematic structure of multiple receiving coils and a 
reception unit in the MRI apparatus according to the third 
embodiment of the present invention. 

Best Mode for Carrying Out the Invention . 
(First embodiment) 

Hereinafter, an MRI apparatus and an MR imaging 
method according to the present invention will be described 
with reference to the drawings. First, a first embodiment of the 
MRI apparatus according to the present invention will be 
described with reference to Fig. 1 to Fig. 7. 

As shown in Fig. 2, the MRI apparatus according to this 
embodiment includes magnet 402 for generating a uniform static 
magnetic field having a predetermined strength to a space 
(hereinafter referred to as "imaging space") having a 
predetermined size for accommodating object 401 to be 
examined, gradient magnetic field coil 403 for generating 
gradient magnetic fields to the imaging space, body coil 1 404 for 
generating a radio-frequency magnetic field to the imaging 
space and detecting an NMR signal generated by object 401, 
and RF probe 405 for detecting the NMR signal generated by 
object 401. Here and hereinafter, "RF" is an abbreviation of 
"Radio Frequency." Although Fig. 2 shows an MRI apparatus of 
horizontal magnetic field type which generates a static magnetic 
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field in a horizontal direction, the present invention can also be 
applied to an MRI apparatus of vertical magnetic field type 
which generates the static magnetic field in the vertical 
direction. 

5 Gradient magnetic field coil 40.3 includes a plurality of 

coils for generating gradient magnetic fields in X-, Y-, and 
Z-directions, each of which generates the gradient magnetic 
field in response to a signal from gradient magnetic field power 
source 409. These gradient magnetic fields are superposed on 

10 the static magnetic field and used for selecting an imaging slice 
and encoding the NMR signal. 

Body coil 404 is designed to generate an RF magnetic 
field (RF pulse) in response to a signal output from RF 
transmission unit 410. Because conductors of body coil 404 

15 are arranged so that its sensitivity distribution becomes 
substantially uniform, it is also used to receive the NMR signal 
used to obtain a reference signal for calculating a sensitivity 
distribution of receiving coil in this embodiment. The RF 
magnetic field is used to excite nuclear spins in the body of 

20 object 404 placed in the imaging space. Meanwhile, the coil 
for obtaining the reference signal for calculating the sensitivity 

distribution of receiving coil may be a dedicated receiving coil 

i 

having a substantially uniform sensitivity distribution, instead of 
the coil for both transmission and reception. 
25 Signals received by body coil 404 and RF probe 405 are 

detected by signal detection unit 406, subjected to a signal 
processing in signal processing unit 407 and then to a 
calculation processing, and converted into image data. The 
image data dealt with in the present invention includes 
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examination image data used in diagnosis on the object and 
sensitivity image data for calculating a sensitivity distribution of 
coil. Among them, the examination image data used in 
diagnosis of the object is displayed as an image on display unit 
5 408. 

_ Gradient magnetic field power source 409, RF 

transmission unit 410, signal processing unit 407, and display 
unit 408 are controlled by control unit 411. A time chart of the 
control of control unit 411 is generally referred to as "pulse 

10 sequence". Object 401 is laid on bed 412, which is designed 
to move between the imaging space inside the magnets and the 
outside of the magnets. 

In the MRI apparatus according to this embodiment, RF 
probe 405 includes two figure-of-eight shaped coils 4051 and 

15 4052 arranged in the inside space of body coil 405, as shown in 
Fig. 3. The two receiving coils 4051 and 4052 are arranged 
opposite at a predetermined interval in the Z direction along an 
X-Y plane, which receive the NMR signal from the object placed 
inside the space surrounded by these two receiving coils 4051 

20 and 4052. Meanwhile, in the figures, elements and 
processings given a four-digit number having an identical first 
three digits represent an identical element and processing, and 
the number of their last one digit is given to distinguish each of 
plural identical elements and processings. Accordingly, they 

25 are hereinafter equivalently treated in this specification. 

Meanwhile, RF probe 405 is constructed with a technique 
referred in the field of MRI as "multiple receiving coils" or 
"phased array coil" including a plurality of receiving coils. Two 
receiving coils 4051 and 4052 J of multiple receiving coils 301 
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are connected with pre-amplifiers 3021 and 3022 as shown in 
Fig. 4. The reception unit for measured signals received by 
multiple receiving coils 301 has signal detection unit 406 and 
signal processing unit 407. In signal detection unit 406, two 
5 AD conversion/quadrature detection circuits 3031 and 3032 are 
arranged in parallel, each of which is connected to outputs of 
pre-amplifiers 3021 and 3022. Signal processing unit 407 is 
designed to generate an MR image from signals received by 
receiving coils 4051 and 4052 by utilizing a Fourier 

10 transformation method, a back projection method, a wavelet 
transformation method, and the like, which includes signal 
combination processing unit 304 for performing the combination 
processing on those image signals. 

The imaging method of thus constructed MRI apparatus 

15 according to this embodiment will be" described. The pulse 
sequence for the imaging may be a general spin-echo pulse 
sequence or a gradient-echo pulse sequence. Furthermore, it 
may be a spin-echo type multi-slice pulse sequence or a 
gradient-echo type multi-slice sequence. For example, in the 

20 gradient-echo pulse sequence shown in Fig. 5, nuclear spins are 
excited by applying thereto RF pulse 601 and slice-selecting 
gradient magnetic field pulse 602. RF pulse 601 applied here 
is provided a frequency on the imaging slice position 
determined depending on a magnetic field gradient of the 

25 slice-selecting gradient magnetic field superposed on the static 
magnetic field, and a frequency band determined depending on 
a slice thickness. After the nuclear spins inside the slice are 
excited, the excited nuclear spins are given positional 
information in a phase direction by the application of phase 
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encoding gradient magnetic field pulse 603. After that, echo 
signal 606 is measured during the application of readout 
gradient magnetic field pulse 604, i.e. during sampling window 
605. This pulse sequence is repeated plural times at time of 
5 repetition (TR) 607. TR 607 is, for example, an interval 
between applications of RF pulse 6011 in a first pulse sequence 
and of RF pulse 6012 in the next pulse sequence. By changing 
an amount of phase encoding gradient magnetic field 603 
(product of strength and application time) at every repetition of 

10 pulse sequence, the phase encoding amount applied to the 
nuclear spins is varied, and echo signals from the nuclear spins 
applied the various phase encoding amounts are detected.. 
This process is repeated as much as the phase encoding step 
number determined in accordance with the pixel number in the 

15 phase encoding direction of the image to be examined, and the 
signals in number determined in accordance with the pixel 
number in the readout direction are sequentially measured in 
sampling window 605. Within image acquisition time 608 a set 
of echo signals necessary for reconstructing one image is 

20 measured. That is, in the conventional imaging method, the 
number of echo signal is determined as 64, 128, 256, 512, or 
the like in consideration of spatial resolution of field of view 
(hereinafter referred to as "FOV") in the phase encoding 
direction. Further, the echo signals are obtained as 

25 time-sequential sampling data usually including 128, 256, 512, 
1024, signals or so in consideration of spatial resolution of FOV 
in the readout direction. Those measured data are stored in a 
two-dimensional memory space (hereinafter referred to as 
"k-space") with the phase direction and the readout direction, 
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read out, and two-dimensionally Fourier-transformed, and one 
MR image is reconstructed. 

According to this embodiment, a plurality of adjoining 
slices of, e.g. m in number are imaged on the examining portion 
5 of object 401. Accordingly, the above pulse sequence is 
sequentially, executed on each of those slice positions. 
Meanwhile, the multi-slice imaging method may be employed 
here as above. 

However, in the MRI apparatus according to this 

10 embodiment, the pulse sequence is modified for executing the 
parallel imaging. That is, according to this embodiment, 
pre-measurement of MR signal is performed with two receiving 
coils, i.e. body coil 404 and receiving probe 405 to rapidly 
perform the parallel imaging. In the MR signal 

15 pre-measurement, first, a pulse sequence of pre-measurement 
in which the size of phase encoding matrix, i.e. the number of 
phase encoding step is reduced to be one of several fractions 
of that in the usual imaging is executed to acquire sensitivity 
image data for n slices included in the examining portion of the 

20 object. Among the sensitivity image data for n slices, the data 
for at least two slices corresponds to the imaging slice position, 
where the slice number is determined as n<m. Further, to 
acquire the sensitivity image data for n slices, the pulse 
sequence shown in Fig. 5 may be applied to the slice positions 

25 where the sensitivity image data is acquired, or the multi-slice 
pulse sequence may be applied so as to measure signals from 
the slice positions where the sensitivity image data is acquired. 

Next, a pulse sequence of imaging (real imaging) in 
which the phase encoding steps in the k space are thinned out 
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is executed to acquire image data for m slices (m>n) of the 
examining portion of object 401 using receiving coils 405t and 
4052 forming multiple receiving coils 405. According to this 
embodiment, the image data for m slices are acquired from 
5 sequentially adjoining slice positions. Meanwhile, to acquire 
the image data for m slices, the multi-slice pulse sequence or a 
pulse sequence for three-dimensional measurement may be 
used instead of the pulse sequence shown in Fig. 5. 

After that, the sensitivity image data acquired by body 

10 coil 404 and receiving probe 405 and the image data acquired 
by receiving coils 4051 and 4052 of receiving probe 405 are 
respectively subjected to two-dimensional Fourier 
transformation and converted into a sensitivity image and an 
examination image. An aliasing artifact is generated in the 

15 examination image reconstructed by the two-dimensional 
Fourier transformation, which is removed with the following 
processings. 

That is, the sensitivity distribution data of receiving coils 
4051 and 4052 are calculated from the sensitivity image 

20 obtained with the pre-measurement pulse sequence, an 
interpolation processing is performed using the calculated 
sensitivity distribution data to calculate the sensitivity 
distribution data corresponding to each slice position of the 
examination image, the sensitivity distribution of receiving coil 

25 is calculated using the calculated sensitivity distribution data, 
and a matrix calculation is performed using this sensitivity 
distribution and the examination image, whereby the aliasing 
artifact on the examination image is removed. 

In the real measurement, the repetition number of pulse 
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sequence is reduced by thinning out the intervals between 
phase encoding steps at a fixed rate. This thinning ratio is 
called "acceleration factor" of the parallel imaging. For 
example, if every other phase encoding steps are thinned out, 
5 the acceleration factor is 2. That is, when the phase encoding 
step number is 64, all data of 1st to 64th phase encodes are 
usually measured. Meanwhile, in the parallel imaging with the 
acceleration factor of 2, data of 1st, 3rd, 5th, .... and 63rd 
phase encodes are measured. At this time, the thinning out 

10 ratio is 1/2. Control unit 411 according to this embodiment 
executes the real-measurement pulse sequence in which the 
phase encoding steps are thus thinned out to acquire the imager 
data for m slices with receiving coils 4051 and 4052. 

Here, the aliasing of image in the parallel imaging will be 

15 described with reference to Fig. 6. Fig. 6(a) shows the case 
where imaging is done with usual phase encoding steps (usual 
imaging), the number of which is not reduced. Signals 2021, 

2022 and 2027 obtained with each phase encoding amount 

are arranged in the k space to form measured data 201 for one 

20 image. By Fourier-transforming it, an image without aliasing is 
obtained as shown in Fig. 6(c). 

On the other hand, as shown in Fig. 6(b), in the parallel 
imaging measurement is performed while the interval between 
the phase encoding steps are e.g. doubled in the same FOV as 

25 in the usual imaging to thin out the phase encoding steps. In 
this case, data are measured on every other lines 2041, 2042, 
and 2043 in a ky direction of the k space, and data on lines 
2051, 2052, .... 2053, and 2054 are not measured. Accordingly, 
the amount of measured data becomes half and the measured 
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FOV also becomes 1/2. An MR signal from a portion of the 
object out of the measured FOV becomes an aliasing 
component. When the k space including measured data 2041, 
2042, and 2043 is Fourier transformed, image 208 in which 
5 an aliasing artifact is generated in the phase encoding direction 
is obtained as shown in Fig. 6(b). That is, object's image 2061 
indicated in upper portion 2071 of image 207 and object's image 
2062 in lower portion 2072 as shown in Fig. 6(c) become 
aliasing and appear as aliasing image 208 shown in Fig. 6(d). 

10 Thus generated aliasing in image 208 can be removed 

with a signal processing method disclosed in, e.g. SENSE: 
Sensitivity Encoding for Fast MRI (Klass P. Pruessmann et. al), 
Magnetic Resonance in Medicine 42: 952-962(1 999), etc. 

Here, a calculation method of aliasing removal will be 

15 described. If the image matrices in the x and y directions are 
X and Y respectively, a pixel value on a coordinate (x, y) (1^x 
^ X, 1^y^Y) in the image is si (x, y). Here, a subscript "i " is 
the number of receiving coil, where 2^i^N. In the case of 
Fig. 6(d), since the phase encoding steps are thinned out to be 

20 half, the matrix in the phase encoding direction of the image 
after subtraction is Y' = Y/2. When a coordinate of the image of 
Fig. 6(d) is (x, y') (1 ^ y' ^ Y'), its pixel value si'(x, y') is 
expressed by formula (1) since two portions 2071 and 2072 of 
original image 207 are superposed on each other. Here, a is a 

25 fixed number. 

s;^y^ s ' (x - y ' )+s ' (x - y ' +r) (D 

a 

Next, the sensitivity distribution of receiving coil and the 
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image will be described. When a two-dimensional sensitivity 
distribution of ith receiving coil is ci(x, y), received signal si(x, 
y) is a product of sensitivity distribution ci(x, y) of receiving coil 
and proton, density distribution p(x, y) of the object, which can 
be expressed by formula (2). 



s,(x,y) = c i (x,y)xp(x,y) 



(2) 
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Here, formula (1) can be expressed by formula (3) using 
formula (2). 



15 



20 



• t/ „ c,(x,/)x A .(x,y) Cl (x,y'+Y')xp(x,y'+Y') 

,(*,/)=— + 

a a 

=-y;c l (x,y+o-i)xr)xp(x,y+a-i)xn (3) 

a % 



When it is given for easy understanding that: 



S, =s\ (*,/) 

C, =c,(x,y+0-l)x7') 

p. =^,y+o--i)xy) 



(4) 



formula (3) is expressed as formula (5). 
S,=-5LC,xPj (5) 
25 Formula (5) can be expressed by formula (6) as a matrix of Nx2. 
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From formula (6), if sensitivity distribution C ti of receiving . coil 
is found, a nucleate factor in measuring the object, e.g. proton 
density distribution Pj can be calculated by calculating the 
inverse matrix. 

Similarly, a general formula in performing imaging at 
M-fold speed using N coils is can be expressed by formula (7), 
where Y' = Y/ M and 1^y'^Y\ Here, character b is a fixed 
number. 
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15 In this manner, by performing the matrix calculation of 

formula (7) on measured data s'i(x, y') using sensitivity 
distribution Cjj of receiving coil, density distribution Pj of 
nuclear spin in the measuring portion is calculated as image 
data in which the aliasing is removed. Meanwhile, a 

20 relationship between the receiving coil number (N) and the 
acceleration factor (M) in the parallel imaging is mathematically 
N ^ M. Further, the sensitivity distribution also may be 
calculated by low-pass filtering the data in a low-frequency 
band of the k space. 

25 In the parallel imaging method of MRI apparatus 

according to this embodiment, sensitivity distribution Cjj of each 
of receiving coils 4051 and 4052 is calculated on the basis of 
sensitivity image data generated from the measured data 
obtained in advance of the low-frequency band in the phase 
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encoding direction of the k space in the pre-measurement pulse 
sequence. 

The sensitivity image data for calculating sensitivity 
distribution Cij of receiving coils 4051 and 4052 are acquired by 
5 executing the pre-measurement sequence. In the 

pre-measurement pulse sequence, the interval between the 
phase encoding steps are not thinned out but the s.ame as in 
the usual imaging. However, by reducing the phase encoding 
step number to be one of several fractions of that in the usual 

10 imaging, the repetition number of pulse sequence is reduced. 
For example, when the phase encoding step number in the k 
space is 64, in the usual imaging all data of 1st to 64th phase 
encode steps are measured. Meanwhile, in the 

pre-measurement sequence of parallel imaging, the phase 
■ 15 encoding step number is set to be the same as in the usual 
imaging, and the data on 16 steps including 25th to 40th phase 
encode steps are measured in one portion of the k space, e.g. 
in the central portion (low-frequency area) of the k space. In 
this case, since the phase encoding step number becomes one 

20 fourth of that in the usual imaging, the time for 
pre-measurement is shortened. Control unit 411 executes the 
pre-measurement pulse sequence to acquire the data to be 
measured of the central portion being about one fourth of the k 
space so as to acquire the sensitivity image data for n slices of 

25 body coil 404 and receiving coils 4051 and 4052. In this 
manner, the sensitivity image data for one slice are measured 
with the small phase encoding step number and the sensitivity 
image data are acquired with the smaller number of slice 
position than that in the examination image, whereby the total 
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imaging time of parallel imaging is shortened. 

Next, a procedure of signal processing performed in the 
signal processing unit of MRI apparatus according to this 
embodiment will be described. 
5 As shown in Fig.1, signal processing unit 407 according 

to this embodiment performs slice interpolation processings 
1011, 1012, and 1013 on sensitivity image 701 of n in number 
of body coil 404 and on sensitivity images 702 and 703, each of 
n in number, respectively of receiving coils 4051 and 4052 of 

10 multiple receiving coils 301 obtained by executing the 
pre-measurement pulse sequence for sensitivity image data 
acquisition, and thus generates an interpolated sensitivity 
image so that all sensitivity images corresponding to the 
respective slice positions of examination image complete. 

15 Next, sensitivity distribution data 707 and 708 respectively of 
receiving coils 4051 and 4052 are calculated by sensitivity 
distribution calculation processings 7061 and 7062 using a 
whole sensitivity image including the sensitivity image 
generated by slice interpolation processing 1011 of body coil 

20 404 and that including the sensitivity image generated by slice 
interpolation processings 1012 and 1013 of receiving coils 4051 
and 4052 of multiple receiving coils 301. Calculated sensitivity 
distribution data 707 and 708 of receiving coils 4051 and 4052 
are stored into a memory (not shown) mounted in control unit 

25 411. 

The above sensitivity distribution calculation processings 
7061 and 7062 are the processing of dividing the sensitivity 
images of receiving coils 4051 and 4052 by the sensitivity 
image of body coil 404 while slice positions are matched. 
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After that, the real-measurement pulse sequence for the 
examination image data acquisition is executed, a determinant 
is formed by matrix forming processing 709 on the basis of 
obtained examination images 704 and 705, each for m slices, of 
5 receiving coils 4051 and 4052 of multiple receiving coils 301 
and sensitivity distribution data 707 and 708, inverse matrix 
calculation processing 710 of the determinant is performed, and 
diagnostic image 711 without an aliasing is resultingly obtained. 
By the way, sensitivity images 701, 702, and 703, each 

10 for n slices, obtained by executing the pre-measurement are 
sensitivity images 7011, 7014, .... and 701m of body coil 404, 
sensitivity images 7021, 7024, and 702m of receiving coil 
4051 of multiple receiving coil 301, and sensitivity images 7031, 
7034, 703m of receiving coil 4052 of multiple receiving coils 

15 301, as shown in Fig. 7. 

Meanwhile, examination images 704 and 705 including an 
aliasing artifact obtained by executing the real-measurement 
pulse sequence are examination images 7041 to 704m of 
receiving coil 4051 of multiple receiving coils 301 and 

20 examination images 7051 to 705m of receiving coil 4052 of 
multiple receiving coils 301. In Fig. 7, ch.1 corresponds to 

receiving coil 4051 of multiple receiving coils 301 and ch.2 

i 

corresponds to receiving coil 4052 of multiple receiving coils 
301. 

25 Each of those sensitivity images 701, 702, and 703 and 

examination images 704 and 705 is given a slice number 
corresponding to its slice position, and images having an 
identical slice number are those on an identical slice position. 
The slice number corresponds to the last one digit of image 
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number. For example, sensitivity images 7011, 7014, and 
701m of body coil 404 are respectively given the slice numbers 
1, 4, and m. Further, slice intervals of the images obtained 
in the pre-measurement sequence and the real-measurement 
5 pulse sequence are evenly determined. 

That is, in the pre-measurement pulse sequence, 
measurement is performed on every three slices as on slice 
numbers 1, 4, and m, and sensitivity image data 701 of body 
coil 404 and sensitivity image data 702 and 703 of receiving 

10 coils 4051 and 4052 are not acquired on slice numbers 2, 3, 5, 
and m-1. Accordingly, in examination images 7042, 7043, 
7045, and 704m-1, and 7052, 7053, 7055, and 705m-1 of 
slice numbers 2, 3, 5, and m-1, the aliasing artifact cannot 
be removed without sensitivity distribution information on those 

15 slice positions. 

Therefore, the signal processing unit according to this 
embodiment calculates sensitivity images 7012, 7013, 7015, 
and 701m-1 by performing slice interpolation processing 1011 
from sensitivity images 701 1, 7014, and 701m of body coil 

20 404, similarly calculates sensitivity image 7022, 7023, 7025, 

and 702m-1 by performing slice interpolation processing 1012 
from obtained sensitivity images 7021, 7024, and 702m of 
receiving coil 4051, and calculates sensitivity images 7032, 
7033, 7035, and 703m-1 by performing slice interpolation 

25 processing 1013 from sensitivity images 7031, 7034, and 
703m of receiving coil 4052. That is, reference number 501 in 
Fig. 7 represents a set of sensitivity image necessary for 
removing the aliasing artifact in the examination image 
including the sensitivity images acquired in the 
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pre-measurement pulse sequence and the sensitivity image data 
acquired by slice interpolation processings 1011, 1012, and 
1013, and reference number 502 represents a set of 
examination image including the aliasing artifact acquired by 
5 the real-measurement pulse sequence. Further, white 
sensitivity images in sensitivity image set 501 represent those 
actually acquired by in the pre-measurement pulse § sequence 
and black ones represent those generated by slice interpolation 
processing 101. 

10 In sensitivity distribution calculation processing 706 

performed in the MRI apparatus according to this embodiment, 
sensitivity distribution ci(x, y) is calculated from image data si(x, 
y) of each receiving coil i on the basis of image data sc(x, y) of 
the body coil, e.g. by relativising it using formula (8): 

15 

C;(x,y) = s i (x,y)xsc(x,y) (8) 

Using thus calculated sensitivity distribution data 707 
and 708 and acquired examination image data 7041 to 704m 

20 and 7051 to 705m respectively of receiving coils 4051 and 4052, 
a determinant of formula (7) is formed by matrix forming 
processing 709. After that, diagnostic image 711 in which an 
aliasing artifact is removed is obtained with reverse matrix 
calculation processing 710. 

25 In a conventional MRI apparatus, although the phase 

encoding step number is set to be one of several fractions of 
that in the usual imaging in the pre-measurement pulse 
sequence, imaging for the sensitivity image data and that for 
the aliasing image data are performed for the same number of 



slices and on the same slice positions. 

However, in the parallel imaging method, the time for 
executing the pre-measurement occupies a large proportion of 
the total imaging time. Especially, in the multi-slice and 
three-dimensional imaging, because the time of repetition (TR) 
and the repetition number of times of pulse sequence increase 
as the number of imaging slice is large, the execution time of 
the pre-measurement pulse sequence is a major matter in 
consideration of the efficiency of examination. 

Meanwhile, in the MRI apparatus according to this 
embodiment, sensitivity image data 701, 702, and 703 are 
acquired for every several slices, e.g. for every three slices to 
reduce the slice number of the imaging sensitivity image data to 
be m being smaller than the slice number m of examination 
image. Accordingly, the repetition number of the 

pre-measurement pulse sequence is reduced, so that the 
imaging time can be shortened for thinned sensitivity image 
data 701 , 702, and 703. Further, when this embodiment is 
applied to three-dimensional measurement, the slice encoding 
step number in the pulse sequence for sensitivity distribution 
acquisition can be reduced. 

Further, in the pre-measurement pulse sequence 
according to this embodiment, since the sensitivity image data 
is acquired while the slice intervals are evenly thinned out, the 
sensitivity image data generated by slice interpolation 
processing 101 are accurately calculated with a simple 
calculation. In addition thereto, in the pre-measurement pulse 
sequence according to this embodiment, setting and control of 
pulse sequence in control unit 411 is simplified since the slice 
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number is thinned out at even intervals to acquire the 
sensitivity image data. 

Further, the sensitivity distributions of receiving coils 
4051 and 4052 are stored in the memory of control unit 411. 
5 Accordingly, when it is necessary to image again 1st to mth 
slices while slightly shifting the slice positions thereof, it is 
possible to calculate the sensitivity distribution of each coil on 
the slice positions to be re-measured by using the sensitivity 
distribution stored in the memory, without re-measuring them. 
10 Therefore, the time for re-measurement can be shortened. 

(Second embodiment) 

Second embodiment of MRI apparatus according to the 

present invention will be described with reference to Fig. 8 and 
15 Fig. 9. Because the same functions as in the first embodiment 

are given the same reference marks, only the difference with 

the first embodiment will be hereinafter described. 

The difference between this embodiment and the first 

embodiment is, as shown in Fig. 8, that slice interpolation 
20 processing 801 is performed after sensitivity distribution 

calculation processing 706 is performed using sensitivity image 

701 of body coil 404 and sensitivity images 702 and 703 of 

I 

receiving coils 7051 and 4052 of multiple receiving coils 301 
acquired in the pre-measurement pulse sequence for the 
25 sensitivity image data acquisition. Therefore, in the signal 
processing unit according to this embodiment, sensitivity 
distribution data 7061 and 7062 of receiving coils 4051 and 
4052 are calculated from sensitivity image 701 of body coil 404 
and sensitivity images 702 and 703 of receiving coils 4051 and 
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4052, the slice interpolation processing 801 is performed next 
using those sensitivity distribution data 706.1 and 7062 to 
acquire sensitivity distribution data 707 and 708 of receiving 
coils 4051 and 4052, matrix forming processing 709 is 
5 performed after that from those sensitivity distribution data 707 
and 708 and examination images 704 and 705, and the aliasing 
artifact in examination images 704 and 705 are removed by 
performing inverse matrix calculation processing 710 on the 
formed matrix . 

10 Here, sensitivity distribution data of receiving coils 4051 

and 4052 calculated by sensitivity distribution calculation 
processing 706 and those generated by slice interpolation 
processing 801 will be described. 

Sensitivity distribution data calculated by sensitivity 

15 distribution calculation processing 706 are, as shown in Fig. 9, 
sensitivity distribution data 7071, 7074, and 707m of 

receiving coil 4051 of multiple receiving coils 301 and 
sensitivity distribution data 7081, 7084, and 708m of 

receiving coil 4052. 

20 In the pre-measurement pulse sequence according to this 

embodiment, the number of measured slices are reduced by 
setting measuring slice positions as every three slices as in the 
first embodiment. That is, the sensitivity image data of slice 
numbers 1, 4, .... and m are acquired and those of slice 

25 numbers 2, 3, 5, and m-1 are not acquired. The aliasing 
artifact appearing in examination images 7042, 7043, 7045, 
and 704m-1 and 7052, 7053, 7055, and 705m-1 cannot be 
removed without the sensitivity distribution data of slice 
numbers 2, 3, 5, and m-1. 



Therefore, in the signal processing unit according to this 
embodiment, sensitivity distribution data 7072, 7073, 7075, 
and 707m-1 are calculated by performing slice interpolation 
processing 8011 on acquired sensitivity distribution data 7071, 
7074, and 707m of receiving coil 4051, and sensitivity 

distribution data 7082, 7083, 7085, and 708m-1 are 

calculated by performing slice interpolation processing 8012 on 

sensitivity distribution data 7081, 7084 and 708m of 

receiving coil 4052. 

In other words, reference number 503 in Fig. 9 represents 
a data set necessary for removing the aliasing artifact in the 
examination image, including sensitivity distribution data 
calculated by sensitivity distribution calculation processings 
7061 and 7062 and those obtained by slice interpolation 
processings 8011 and 8012. The white sensitivity image in 
data set 503 represents the sensitivity distribution data 
calculated by sensitivity distribution calculation processings 
7061 and 7062, and the black sensitivity distribution data 
represent sensitivity distribution data generated with by 
interpolation processing 801. 

In this manner, signal processing unit 407 performs slice 
interpolation processings 8011 and 8012 after performing 
sensitivity distribution calculation processings 7061 and 7062 
on the basis of sensitivity image data 701, 702, and 703 
acquired in the pre-measurement pulse sequence, forms 
sensitivity distribution data 7071 to 707m and 7081 to 708m 
corresponding to each slice position of examination images 
7041 to 704m and 7051 to 705m, and after that removes the 
aliasing artifact in the image by matrix forming processing 709 
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and inverse matrix calculation processing 710. 

Thus, according to this embodiment, since slice 
interpolation processing 801 is performed after calculating 
sensitivity distributions of receiving coils 4051 and 4052 on the 
basis of sensitivity images 701, 702, and 703 acquired in the 
first pulse sequence, it is unnecessary to perform the slice 
interpolation processing on sensitivity image 701 of body coil 
404. Therefore, in comparison with First embodiment where 
slice interpolation processing 1011 is performed on sensitivity 
image 701 of body coil 404, it is possible to simplify sensitivity 
distribution calculation processing 706 and shorten the time for 
image reconstruction. 

(Third embodiment) 

Next, the third embodiment of MRI apparatus according 
to the present invention will be described with reference to 
Fig. 10 and Fig. 11. Meanwhile, in this embodiment, too, the 
same function as in First embodiment is given the same 
reference mark and the description thereof will be omitted. 
The structure different from First embodiment and the 
characteristic part will be described hereinafter. 

The difference between this embodiment and First 
embodiment lies in that the multiple receiving coils include four 
receiving coils. That is, as shown in Fig. 10, multiple receiving 
coils 301A forming an RF probe in the MRI apparatus according 
to this embodiment include four figure-of-eight shaped receiving 
coils 4051 to 4054, which are arranged opposite on the X-Y 
plane and on the T-Z plane at a predetermined distance. 

Further, as shown in Fig.11, four receiving coils 4051 to 
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4054 of multiple RF receiving coils 3.01 A according to this 
embodiment are respectively connected to preamplifiers 3021 to 
3024. In signal detecting unit 406A, four AD 

conversion/quadrature detection circuits 3031 to 3034 are 
5 arranged in parallel, which are respectively connected to 
outputs of preamplifiers 3021 to 3024. 

By the way, when the parallel imaging is executed, a 
portion of low sensitivity is generated in some cases depending 
on the arranging relationship and the sensitivity distribution of 
10 the plural receiving coils. In this case, the matrix calculation 
of formula (7) may diverge to generate an artifact or deteriorate 
image quality. 

In the MRI apparatus according to this embodiment, 
among four receiving coils two optimum combinations of 

15 receiving coils is formed on the basis of an imaging cross 
section and the phase encoding direction. Here, "optimum 
combination of receiving coil" is a combination with which the 
matrix calculation does not diverge. For example, there may 
be such conditions that there is not a portion of low sensitivity 

20 distribution in the phase encoding direction, that sensitivity 
distributions of receiving coil groups are not identical, and the 
like. 

As shown in Fig. 10, two receiving coils 4051 and 4052 
according to this embodiment are arranged opposite in the Z 
25 direction on the X-Y plane at a predetermined distance. 
Further, two receiving coils 4053 and 4054 are oppositely 
arranged in the X direction on a plane parallel to the Z-Y plane 
at a predetermined distance. Among those four receiving coils, 
receiving coils 4051 and 4053 are combined to form receiving 
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coil group 4055, and receiving coils 4052 and 4054 are 
combined to form receiving coil group 4056. When the 
receiving coils are thus arranged and combined, the combined 
sensitivity distribution data of receiving coils 4051 and 4053 
5 and that of receiving coils 4052 and 4054 become different from 
each other and the number of low sensitivity portions decreases, 
whereby the matrix calculation does not easily diverge. 

Because this optimum coil combination is determined 
depending on the imaging cross section and the phase encoding 
10 direction if the coil arrangement is fixed,, the apparatus is 
desirably constructed so that the optimum combination is 
automiatically selected when those imaging conditions are 
determined. 

In signal processing unit 407A, two sets of sensitivity 
15 image data measured by the two selected combinations of 
receiving coils are picked up, and those sensitivity image data 
of each set are combined. Similarly, two sets of examination 
image data measured by the two selected combinations of 
receiving coils are picked up and those examination image data 
20 of each set are combined. The sensitivity, image data of 
receiving coil group 4055 generated by thus combining those 

data correspond to sensitivity image data 7021, 7024, and 

702m in ch.1 of Fig. 7, and those of receiving coil group 4056 
correspond to sensitivity image data 7031, 7034, .... and 703m 
25 in ch.2. 

Similarly, the combined examination image data of 
receiving coil group 4055 correspond to examination image data 
7041 to 704m in ch.1 of Fig. 7 and those of receiving coil group 
4056 correspond to examination image data 7051 to 705m in 
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ch.2. 

That is, the sensitivity image data acquired by receiving 
coils 4051 to 4054 are combined in signal combination 
processing unit 304A into sensitivity image data 702 and 703 in 
5 ch.1 and ch.2 corresponding to the sensitivity distribution of 
each predetermined receiving coil group, and also made into 
reference data of the matrix calculation for removing the 
aliasing artifact in examination image data 704 and 705 of ch.1 
and ch.2 combined in each receiving coil group. 

10 Here, because the sensitivity image data and 

examination image data acquired by receiving coils 4051 to 
4054 are complex numbers, the combination processing based 
on , e.g. complex addition is done in signal combination 
processing unit 304A. Alternatively, however, the combination 

15 processing may be done by calculating a root sum square of two 
sets of sensitivity image data and that of two sets of 
examination image data. The latter case brings an advantage 
of higher SN ratio in comparison with the complex addition. 
Meanwhile, in the procedure of signal processing shown in 

20 Fig . 1 1 , the signals are combined after the AD conversion. 
However, combining processing may be done on analog signals 
before being subjected to the AD conversion, and the AD 
conversion is performed afterward on them. 

In this manner, in the rapid parallel imaging using the 

25 multiple receiving coils in the MRI apparatus according to this 
embodiment, an optimum combination of receiving coils is 
selected considering the imaging cross section and the 
encoding axis, and the matrix calculation is performed on the 
combined measured data of the receiving coil combination, 
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whereby image deterioration caused by the matrix calculation is 
avoided and a good image is obtainable. Matrix [GxM] is 
calculated in this matrix calculation, where the thinning out 
ratio of signal is 1/M and the combination number of receiving 
5 coils is G (here, G = 2, N, which corresponds to N in formula 
(7)). According to Third embodiment, since four receiving coils 
4051 to 4054 are combined into two receiving coil groups (G=2) 
and the matrix calculation is done on the combined measured 
data of each of two receiving coil groups, the matrix is [2x2], 

10 and thus the calculation can be simplified. 

Further, even when the imaging cross section and the 
encoding axis are changed, a stable image without an artifact 
and deterioration is obtainable. Moreover, it is possible to 
reduce a burden of the matrix calculation and rapidly 

15 reconstruct an image. 

In this embodiment, when the first combination and the 
second combination of receiving coils are selected, one 
receiving coil may be overlapped in both the first and the 
second combinations. That is, the plural combinations of 

20 receiving coils may include a common receiving coil. Since the 
plural combinations of receiving coils may thus share a part of 
receiving coils, the coil combinations with which an optimum 
sensitivity distribution is obtainable can be more flexibly 
selected and a better quality image is obtainable. 

25 The present invention is not limited to the above 

disclosed embodiments and may employ various embodiments 
in consideration of the gist of the invention. For example, 
although First to Third embodiments are described with regard 
to the case where the pre-measurement pulse sequence for the 
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sensitivity distribution acquisition is performed on every three 
slices, it may be done on every arbitrary number of slices. 

Further, although First to Third embodiments are 
described with regard to the case where the slice intervals in 
5 acquiring the image data are equivalent, they may be changed 
to acquire the image data in the MRI apparatus, according to the 
present invention. 

Further, although two receiving coils are used according 
to First and Second embodiments and four receiving coils are 

10 used according to Third embodiment, the number of receiving 
coils is not limited to two and four. That is, it may be eight or 
other number. However, to simplify the circuitry of the 
apparatus and to reduce the object's uncomfortable feeling, it is 
considered that the number of receiving coils is desirably 

15 smaller. 

Further, although First to Third embodiments are 
described with regard to the case where the interval between 
the phase encoding steps in the pulse sequence for calculating 
the sensitivity of receiving coil is twice as in the usual imaging, 

20 or the thinning ratio is 1/2, the step interval may be set as three 
times, four times, or more. In this case, it is necessary to set 
number G of selected combination of receiving coils to be the 
same number as interval M of the phase encoding steps or more 
However, G = M is preferable for simplifying the matrix 

25 calculation. 

Further, although First to Third embodiments are 
described with regard to the case of employing the gradient 
echo sequence, the type of the parallel imaging and the 
measurement method are not limited. For example, it is 
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possible to employ the SE pulse sequence, the FSE pulse 
sequence, the EPI pulse sequence, the SSFP method, and the 
spiral measurement method. 

Further, the MRI apparatus according to the present 
5 invention is applicable not only to the multi-slice measurement 
method of these embodiments but also to the three-dimensional 
measurement method. In this case, the data may be thinned 
out not in the phase encoding direction but in the slice encoding 
direction, or in both the phase encoding direction and the slice 

10 encoding direction in combination. Alternatively, the 
combination of receiving coil may be selected so that the 
combination of sensitivity distributions is optimum in the axis 
direction in which the encoding steps are thinned out. 

Further, in the signal processing done in the signal 

15 processing unit according to the present invention, slice 
interpolation processings 101 and 801 may be, e.g. a spline 
interpolation processing, an interpolation processing based on a 
Sine function, a function-based fitting processing, and the like. 

Further, in the MRI apparatus according to the present 

20 invention, an approximate sensitivity distribution of receiving 
coil can be calculated by generating a sensitivity image data 
having a uniform sensitivity distribution by combining sensitivity 

i 

image data acquired by plural receiving coils instead of the 
body coil according to those embodiment, and dividing the 
25 image data acquired by each receiving coil by thus combined 
sensitivity image data. However, the sensitivity distribution of 
receiving coil can be accurately calculated by dividing the 
image data acquired by each receiving coil by the image data of 
body coil having a relatively uniform sensitivity distribution 
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described in the above embodiments and calculating the 
approximate sensitivity distribution of receiving coil. Therefore, 
the latter is preferable. 

Further, the multiple receiving coils of the MRI apparatus 
5 according to the present invention are not limited to those 
according to First to Third embodiments, and various multiple 
receiving coils may be utilized, such as ones mentioned in 
"Array Head Coil for Improved Functional MRI" (Christoph 
Leussler),1996 ISMRM abstract p. 249, ones mentioned in 

10 "Helmet and Cylindrical Shaped CP Array Coils for Brain 
Imaging: A Comparison of Signal — to — Noise Characteristics" 
(H. A. Stark, E. M. Haacke), 1996 ISMRM abstract P. 1 41 2 as a 
head QD multiple RF coils for horizontal magnetic field, and 
ones mentioned in "Four Channel Wrap-Around Coil with 

15 Inductive Decoupler for 1.5T Body lmaging"(T. Takahashi et. al), 
1995 ISMRM abstract P. 1418 as an abdominal QD multiple RF 
coils for horizontal magnetic field. The multiple receiving coils 
are designed to enlarge an FOV while maintaining high 
sensitivity of receiving coil to improve the whole sensitivity by 

20 combining signals acquired by each of arranged small-sized 
receiving coils having a relatively high sensitivity. Accordingly, 
it is possible to enlarge the FOV and improve the whole 
sensitivity by combining the reconstructed images of each of 
those receiving coils into one image. 

25 Further, although First to Third embodiments are 

described for the case of utilizing a horizontal magnetic field 
type MRI apparatus in which the static magnetic field direction 
(Z direction) is horizontal, the present invention is not limited 
thereto and also applicable to an MRI apparatus in which the 
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static magnetic field direction is vertical. 

As described above, according to the present invention, 
because the sensitivity distribution of receiving coil for 
removing an aliasing artifact generated in an examination image 
5 obtained in the parallel imaging can be acquired by measuring 
data on fewer slices than imaging slices, it can be acquired 
within a short time. Therefore, the total imaging time in the 
parallel imaging can be shortened. Further, since an 
unacquired sensitivity distribution of the receiving coil on 
10 imaging slices is calculated with the interpolation calculation, 
the calculation can be easily performed. 



